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We explore the phase diagram of two flavour QCD at vanishing chemical potential using dy- 

r — 

namical O (a) -improved Wilson quarks. In the approach to the chiral limit we use lattices with 
^sO a temporal extent of N t = 16 and spatial extent L = 32,48 and 64 to enable the extrapolation to 

the thermodynamic limit with small discretisation effects. In addition to an update on the scans at 
constant fc, reported earlier, we present first results from scans along lines of constant physics at a 
pion mass of 290 MeV. We probe the transition using the Polyakov loop and the chiral condensate, 
as well as spectroscopic observables such as screening masses. 
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1. Introduction 

The unresolved question about the order of the phase transition connected to chiral symmetry 
restoration in the chiral limit of two-flavour QCD is the remaining qualitative issue concerning 
the phase diagram in the {m u d,m s ,T} parameter space at zero chemical potential (see [1] for a 
review). There are two possible scenarios [2, 3]: In the first scenario the chiral critical line reaches 
the m U( j = 0-axis at some tri-critical point mj nc and the transition at Nf = 2 with m ut j = is of second 
order. Then the restoration of chiral symmetry belongs to the SU (2) x SU (2) ~ 0(4) universality 
class [4]. In the second scenario the chiral critical line never reaches the m uc [ = axis and the 
transition remains first order for all values of the strange quark mass. In [2, 3] it was shown that the 
realisation of one of the two scenarios can be linked to the strength of the anomalous breaking of 
the ?7A(l)-symmetry at the transition point in the chiral limit. The strength of the anomaly can be 
probed by looking at correlation functions in scalar and pseudo-scalar channels and the associated 
screening masses. Despite a number of recent studies aiming to extract information about the order 
of the transition in the chiral limit [5, 6, 7, 8] no study has sufficient control over the systematic 
effects. 

In these proceedings we present the current status of our study on the topic, first reported 
in [9, 10], using non-perturbatively G (a) -improved Wilson fermions at Nf = 2. We use lattices 
with a temporal extent of N t = 16 throughout to suppress discretisation effects, in particular the 
effect of the explicit breaking of chiral symmetry introduced by the Wilson term. We aim at sim- 
ulations at several quark masses corresponding to zero-temperature pions with m% < 300 MeV 
with three different volumes for each simulation points. This will eventually enable us to perform 
a scaling analysis with control over the main systematic effects. Here we present results on the 
transition temperatures, screening masses and on the strength of the anomalous breaking of the 
t/A(l)-symmetry. 

2. Setup 

The simulations are done using non-perturbatively ^(a)-improved Wilson fermions [11] with 
two degenerate dynamical quarks and the configurations are generated using DD-HMC [12, 13] 
and MP-HMC [14] algorithms. For more details on the simulation setup see [9]. 

To extract information about the order of the transition in the chiral limit by means of a scaling 
analysis, it is important to control and disentangle the different systematic effects that might distort 
the scaling properties of the results at finite lattice spacing and volume. Of particular importance 
for Wilson fermions in this context is the suppression of the explicit breaking of chiral symmetry 
which shows up as a discretisation effect at finite lattice spacing. Therefore our simulations are 
done on large lattices of the size of 16 x 32 3 , 16 x 48 3 and 16 x 64 3 . The three different volumes 
at each simulation point allow for an extrapolation to infinite volume, and at this large temporal 
extent one can expect discretisation effects to be small (see also [15]). 

We scan in the temperature by varying the bare coupling /3 while N t = 16 remains fixed either 
at fixed hopping parameter K* (for heavier quarks) or along lines of constant physics (for quarks with 
an associated m % < 290 MeV). We set the scale by using an interpolation of the zero-temperature 
results for the Sommer parameter r^/a in the chiral limit obtained within the CLS effort in the 
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Figure 1: Simulation points in the {m M /, r}-parameter space. The black dashed line results from naive 
0(4)-scaling for the two transition temperatures indicated by the grey squares and is only shown to give a 
first impression on the possible location of the critical line. 



region of 5.20 < j8 < 5.50 [16, 17] and its continuum result ro = 0.503 (10) fm obtained in [17]. 
Similarly, the tuning of the bare parameters to lines of constant renormalised quark masses is done 
using the input from results obtained within CLS [17, 18, 19]. 

We extract the critical temperature for the deconfinement transition and the chiral symme- 
try restoration using the real part of the APE-smeared Polyakov loop L$m = Re K^sm)] and the 
subtracted chiral condensate [20, 21] 

< W> sub = 2^- mpcAC j d A x <P(0) P(x)) , (2. 1) 

respectively. Here P(x) is the pseudoscalar density and mpcAC the PCAC mass. The transition 
temperatures are defined by the position of the peak of the associated susceptibilities 

x (0)^N^(0 2 )-(0) 2 ), (2.2) 

where O is any of the observables above. For the time being all those observables are unrenor- 
malised. The error analysis has been done using the bootstrap method with 1000 bins. 



3. Current status 

The current set of simulation points in the {m„j, r}-parameter space 1 is shown in figure 1. 
The set of points consists of two different temperature scans, one at constant K with two volumes 
(red and magenta circles) and another one at a line of constant quark mass of m U( j = 14.5 MeV 
(orange diamonds). The associated parameters are listed in table 1. 

'Here m uc ] is the renormalised quark mass in the MS-scheme at a renormalisation scale of fl = 2 GeV. 
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scan 


Lattice 


DD 


K 


m ud [MeV] 


T [MeV] 


X Up [MDU] 


MDUs 


B1 K 


16 x 32 3 


8 4 


0.136500 




190-275 


~ 10 


~ 20000 


B3 K 


16 x 64 3 


8 x4 3 


0.136500 




240-250 


~ 46 


~ 16000 


CI 


16 x 32 3 


8 4 




14.5 (45) 


175-250 


~ 14 


~ 12000 



Table 1: Scans at N t = 16 at constant K, the ones with subscript K", and at constant renormalised quark 
mass, CI. Listed are the size of the DD-HMC blocks, DD, the temperature range in MeV, the integrated 
autocorrelation time of the plaquette Xjj p and the number of molecular dynamics units, MDUs, used in the 
analysis. The measurements have been done each 4 MDUs. 



The results for scan B1 K have in part already been reported in [9, 10]. The main results 
remain unchanged even with increased statistics and a larger number of simulation points. The 
transition temperature in physical units, listed in table 2, changes slightly due to the updated scale 
determination reported in [17]. Note that the transition temperature is extracted from the peak 
position of a Gaussian fit to the Polyakov loop or condensate susceptibility peak. The uncertainties 
are estimated conservatively by the full spread of points included in the fit. The value given in 
table 2 for scan B1 K is extracted from the Polyakov loop, but the condensate also shows a (weak) 
peak in its susceptibility at similar temperature. The additional simulations points from scan B3 K 
probe the same transition only with twice the volume. The results are in good agreement with the 
results for the smaller volume but the limited number of three simulation points does not allow for 
a reliable determination of the critical temperature. The scan serves as a benchmark for the future 
simulations on larger volumes. 

Scan CI is our first scan along a line of constant quark mass of 14.5 MeV, which corresponds to 
a zero-temperature pion mass of about 290 MeV, and provides first results in the regime relevant for 
the future scaling analysis. The results for the smeared Polyakov loop, the subtracted condensate 
and their susceptibilities are shown in figure 2. As can be seen from the plot, both susceptibilities 
exhibit a peak at a similar temperature. The result given in table 1 is the one extracted from the 
condensate, which agrees with the one extracted from the Polyakov loop within errors. 



scan 


T c [MeV] 


m u d,c [MeV] 


m^c [MeV] 


B1 K and B3 ff 


245 (7)(6) 


45 ( 2)( 2) 


512(30)(15) 


CI 


2H(5)(3) 


14.5 (20)(13) 


287 (25)(13) 



Table 2: Estimates for the transition temperatures, quark and pion masses. The first error bar reflects the 
uncertainty due to the extraction of the peak position. The second error denotes the uncertainty due to 
scale setting and renormalisation. For scan Bl^ and B3k the result for the transition temperature has been 
extracted from the Polyakov loop susceptibility, for scan C2 the value from the subtracted condensate has 
been used. The conversion from quark to pion masses is done using continuum ^PT to NNLO [22] with the 
low energy constants from [19, 23]. 
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Figure 2: Results for the smeared Polyakov loop LgM (left) and the subtracted chiral condensate (YV) S ub 
(right) and their susceptibilities for scan CI. The coloured areas are the transition regions extracted from the 
susceptibility. The curves are the Gaussian fits used to define the transition point. 



4. Screening masses 

Further information about the chiral symmetry restoration pattern in a given scan can be ex- 
tracted from the behaviour of screening masses [24]. In particular, the degeneracy of pseudoscalar 
and scalar screening masses signals the restoration of the anomalously broken L^(l)-symmetry. 
In this section we focus on scan CI and study screening masses in the pseudoscalar (P), scalar 
(S), vector (V) and axial vector (A) channels, measured on the stored configurations (separated by 
40 MDUs) with a point source. 

Figure 3 illustrates the temperature dependence of the screening masses. The *-axis is nor- 
malised to the critical temperature listed in table 2. At 0.84 Tq the screening masses are mostly in 
agreement with the expected splitting patterns of the zero-temperature meson masses. The screen- 
ing mass in the pseudoscalar channel, Mp, starts from a value of Mp/(2n T) = 0.33 (2), which 
means that at this point Mp is a factor of 1 .27 larger than the zero-temperature pion mass in the 
scan. Around Tq it then starts to rise and at 1.16 Tq it is roughly 30 % smaller than the asymptotic 
value 2k T . This is in the ballpark of what has been found for the Wilson action on pure gauge 
configurations [25], but larger than typical results for staggered fermions [26, 27]. Note however, 
that regardless of the discretisation finite volume (here N s /N t = 2) and quark mass effects might 
still give a sizeable contribution around Tq- The V and A-channels are accidentally already close to 
the asymptotic value of 271 T below 7c and mainly fluctuate around this value in the whole interval. 
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Figure 3: Temperature dependence of the screening masses in P, S, V and A channels. The dashed line 
corresponds to the asymptotic value of M = 2k T. 

My and Ma fluctuate independently below Tc and become degenerate above Tc, consistent with 
chiral symmetry restoration. 

The strength of the anomalous breaking of the chiral £/A(l)-symmetry can be assessed via the 
splitting between the screening masses in P and S channel. Below and at Tc the screening mass 
in the 5-channel and shows large fluctuations. Above Tc the signal becomes more stable and Ms 
moves closer to Mp, signalling a weakening of the breaking of £/a(1)- At 1.16 7c the symmetry is 
almost restored. This is in qualitative agreement with the findings from [28] where the symmetry 
has been found to be restored at about 1.25 Tc- 

5. Conclusions 

This proceedings article contains a summary of the status of our ongoing study of the QCD 
deconfinement transition in the chiral limit at Nf = 2. To date there are two scans available with 
N t = 16 and pion masses at the critical point of about 510 and 290 MeV. The critical temperatures 
are 245 and 211 MeV respectively. The scan at the pion mass of 290 MeV is the first scan in the 
region important for a future scaling analysis and is done along a line of constant renormalised 
quark mass of 14.5 MeV. To study the pattern of chiral symmetry restoration, we calculate mesonic 
screening masses. The pseudoscalar screening mass rises from a value close to the zero-temperature 
pion mass at 0.847c towards the asymptotic value of 2nT. At 1.16 7c it differs by roughly 30 % 
from this limit. This is in agreement with the results for screening masses extracted from pure 
gauge theory with Wilson fermions [25], but larger than typical results found in simulations with 
staggered fermions [26, 27]. The results show the expected degeneracy for vector and axial vector 
channels around and above Tc, signalling chiral symmetry restoration. At the same time, the Ua{X)- 
symmetry still appears to be broken. The large lattices used in our study also offer the possibility 
to study plasma properties in terms of temporal correlation functions in the vector channel and the 
associated spectral function (see [29]). A first study has already been reported at conferences and 
we refer to our future publication for the details. 
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